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SECTION  I 


INTRODUCTION 

The  analysis  of  multi conductor  transmission  line  networks  plays 
an  important  role  in  studying  the  electromagnetic  pulse  (EMP)  response 
of  aircraft  or  other  aerospace  systems.  Transmission  lines  may  easily 
pick  up  the  EMP  transient  fields  which  penetrate  into  the  system  through 
slits,  holes,  or  other  apertures  in  its  exterior  skin.  Many  of  these 
lines  are  connected  to  vital  electronics  which  control  the  functioning 
of  the  in-flight  vehicles. 

Recently  Tesche  and  Liu  [Ref.  II  have  developed  a  computer  code, 
denoted  QV7T,  to  analyze  a  general  lossless  multi  conductor  network.  The 
code  does  the  analysis  in  the  frequency  domain.  The  Fast  Fourier  Trans¬ 
form  (FFT)  algorithm  is  used  to  obtain  the  transient  responses.  In  con¬ 
trast  with  direct  time  domain  analysis,  tne  QV7T  code  solves  all  the 
propagation  modes  simultaneously,  by  making  use  of  the  BLT  equations 
[Ref.  21  and  matrix  methods. 

The  purpose  of  this  report  is  to  compare  the  transient  responses 
calculated  from  the  QV7T  code  with  measured  data  for  selected  multi- 
conductor  networks.  The  transient  measurements  include  two  types,  i.e., 
the  laboratory  measurement  on  several  experimental  models  and  the  EMP 
simulation  measurements  of  the  cables  inside  an  actual  aircraft. 

The  laboratory  measurements  have  been  provided  by  Mission  Research 
Corporation,  using  test  facilities  at  the  Air  Force  Weapons  Laboratory 
(AFWL)  [Ref.  3,  4].  Three  different  types  of  networks  will  be  considered 
in  this  report:  a  3-wire  transmission  line  section,  a  T-network,  and  an 
H-network.  The  pulse  responses  and  unit  step  responses  were  included  in 
the  laboratory  measurements. 

The  EMP  simulation  measurement  involves  a  power  cable  tested  in 
the  TACAMO  Add-On  Testing  program  [Ref.  51.  A  seven  wire  power  cable 
from  engine  No.  1  runs  parallel  to  the  edge  of  the  wing  and  picks  up 
the  EMP  energy  from  outside  the  pressure  break.  This  energy  is  carried 
into  the  interior  of  the  aircraft  by  the  power  cable.  Measured  data 
enable  the  Thdvenin  equivalent  circuit  to  be  used  to  represent  the 
sources.  The  short  circuit  currents  of  the  cable  at  the  distribution 
box  were  measured  and  are  compared  with  the  calculated  results. 


SECTION  II 


DATA  COMPARISONS  WITH  LABORATORY  MEASUREMENTS 

From  discussions  with  the  personnel  of  Mission  Research  Corpora¬ 
tion,  we  have  obtained  laboratory  measurements  from  some  simplified 
models  of  multi conductor  transmission  line  networks  [Ref.  3] .  Three 
different  networks  were  provided:  the  3-wire  transmission  line  section, 
the  T-network,  and  the  H-network.  The  wires  of  the  3-wire  section  and 
the  T-network  were  insultated  with  polyethylene,  neoprene  and  rubber. 

The  H-network  was  a  single  wire  network  (over  a  ground  plane)  and 
consisted  of  wires  without  dielectric  cladding.  For  all  the  cases  con¬ 
sidered  here,  the  wires  of  the  network  were  supported  with  styrofoam 
blocks  above  an  aluminum  ground  plane.  Time  domain  ref 1 ectometry  (TDR) 
measurements  were  used  to  determine  the  transmission  line  parameters, 
such  as  the  per-unit-length  capacitance,  and  inductance  matrices. 

The  driving  waveforms  for  these  cases  were  a  pulse  excitation 
for  the  3-wire  line  section  and  the  T-network,  and  a  unit-step  function 
for  the  T-network  and  the  H-network.  The  pulsed  driving  waveform  was 
digitized  and  linearly  interpolated  for  use  in  the  FFT  routines  in  the 
QV7T  program.  The  QV7T  program  has  also  been  modified  to  employ  the 
numerical  driving  waveform  data  of  the  excitation  function  for  each 
wire  of  the  network.  Table  la  shows  the  original  UPDATE  correction 
set,  and  Table  lb  shows  the  additional  UPDATE  set  of  the  QV7T  code 
for  the  interpolation  of  the  driving  waveform. 

The  transient  response  for  each  network  configuration  has  been 
obtained  both  by  calculation  and  by  the  laboratory  measurement.  The 
case  by  case  comparison  of  the  calculated  and  the  experimental  results 
will  be  considered  in  the  following  sections.  Through  the  comparisons 
we  shall  be  able  to  establish  the  confidence  of  the  accuracy  of  the 
computational  techniques,  as  well  as  the  adequacy  of  the  measurement 
procedures. 


TABLE  la.  THE  ORIGINAL  UPDATE  CORRECTION  SET  FOR  QV7T  CODE 


*ID  JUN8 
*D  MAIN. 3 

1  TAPE8 ,TAPE9 ,TAPE10 ,TAPE21 
*D  TYMDRI V . 9 .TYMDRIV . 10 
*1  TYMDRIV. 13 

DIMENSION  TT(20) ,VTT(20) ,C0EF(2.20) 

NAMEL 1ST /WV  FORM/ 1 STE  P ,T0 
NAMELIST /FFTDR/TMAX ,M 
NAMEL 1ST /WVMATCH/NTT , SCAL , TT , VTT 
‘D  DTUBE.14 
‘D  D TUBE. 17 
*1  DTUBE.18 

NAMELIST/LINE PRM/XK,SL,SLEN,IFLAG 
NAMEL IST/SOURCES/XLPOS ,VSRC ,CSRC 
*D  TYMDRIV. 61 .TYMDRIV. 63 
ALPHA1=4.0E6 
ALPHA2=4.76E8 

TM=ALOG ( AL  PHA1 /AL  PHA2 ) / ( AL  PHA1 -AL  PHA2 ) 

A1=1.0/(EXP( -ALPHA1*TM) -EXP( -ALPHA2*TM) ) 

*0  DJUNC.105 
*1  DJUNC.103 

ZL(KR,KC)=0.0 
*1  TYMDRIV. 161 

WRITE  (2)  NLOAD 
*1  TYMDRIV. 195 

WRITE  (2)  NL.NN 
PRINT  100, NL 

100  FORMAT ( 1H0 , 10X , *PR I NT  OUT  OF  LOAD  QUANTITIES  FOR  L0AD*13// 

1 16X , *T IME ( NS ) * , 10X , ‘VOLTAGE ( V ) * , 10X , ‘CURRENT  ( AMPS ) */ ) 

DO  101  1=1, NN 
VV=REAL(VT( I) ) 

V I =REAL ( IT ( I ) ) 

WRITE  (2)  XT ( I ) ,VV,VI 
PRINT  102,XT( I) ,VV ,VI 

102  F0RMAT(1H  ,8X,F12.3,9X,E13.3,9X,E13.3) 

101  CONTINUE 
‘TYMDRIV. 131 

WRITE  (2)  NT,NN1,(VT(I),XT(I),I=1,NN1) 

*D  TYMDRIV. 202, TYMDRIV. 203 

22  F0RMAT(*1  DRIVING  WAVEFORM  FOR  TUBE  *13*  IS  A  DOUBLE  EXPONENTIAL 
1  WITH  THE  FOLLOWING  PARAMETERS*/10X,*A1=E12.3*ALPHA1=*E12.3/  1 

‘D  DJUNC.89 

C..  JZL=1  ..  ZL  IS  A  MATRIX 
*D  JNSCM.156 
*1  JNSCM.150 
NPP=NWW 
‘D  MODAL. 23 

C  DIAGONISTIC  PRINT  OUT  OCCURS  IF  IDBUG.GE.3 
*D  DTUBE.21 

C  LOOP  OVER  EACH  TUBE 


TABLE  lb.  THE  ADDITIONAL  UPDATE  CORRECTION  SET  FOR  INTERPOLATING 

THE  PULSED  DRIVING  WAVEFORM 

*1  TYMDRIV.74 

C  ISTEP  =  3  FOR  ACTURAL  WAVEFORM  EXCITATION 
*1  TYMDRIV.76 

I F( ISTEP . EQ . 3 )  Go  To  301 
*1  TYMDRIV.125 

301  CONTINUE 
C 

C  ACTUAL  WAVEFORM  EXCITATION 
C 

READ  (5.WVMATCH) 

DO  306  1=1, NTT 
VTT ( I ) = VTT ( I ) *SCAL 
306  CONTINUE 

PRINT  304, NT,NTT,(III1,TT(III1),VTT(III1), 1111=1, NTT) 

304  FORMAT ( 1H 1 , 5X ,*DR I V I NG  WAVEFORM  FOR  TUBE  *,13, 

1  *  IS  AN  ACTUAL  WAVEFORM  WITH*, 13,*  SAMPLE  POINTS  * 

2  *AS  F0LL0WS*/5X ,*POINT  N0.*,3X,*TIME(NS)*,5X, 

3  *AMPL ITUDE( VOLTS )*/20(8X , 13 ,6X ,F8 . 2 ,5X ,E13 . 5/ ) ) 

NTT1=NTT-1 

DO  302  I=1,NTT1 

C0EF(1,I)-(VTT(I+1)-VTT(I))/(TT(I+1)-TT(I)) 

C0EF(2,I)=VTT( I)-C0EF(1,I)*TT(I) 

302  CONTINUE 
TTTMAX=TT(NTT) 

11=1 

DO  303  1=1, NN1 
T=DT*FL0AT(I-1)*1.E9 
111=11 

DO  305  II2=II1,NTT1 
IF(T.GT.TT( I 1+1 ) )  11=11+1 
IF(Il.GE.NTT)  11=11-1 

305  CONTINUE 

IF(T.LE.TTTMAX)  VT( I ) =COEF ( 1 , I 1 )*T+COEF( 2 , I 1 ) 

IF(T.GT.TTTMAX)  VT(I)=0.0 
VT(NN1+I)=-VT(I) 

303  CONTINUE 
GO  TO  13 
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A.  3-wire  Transmission  Line  Section 

The  first  case  to  be  discussed  is  a  3-wire  transmission  line 
section.  Figure  1  shows  the  configuration  of  this  line.  Each  wire  is 
terminated  by  a  50ft  resistor  to  ground.  The  pulse  generator  has  an 
effective  resistance  of  25ft,  due  to  the  parallel  combination  of  the 
50ft  load  impedance  and  the  internal  impedance  of  the  generator. 

The  network  topology  for  this  case  is  rather  simple.  There 
are  one  branch  and  two  junctions  at  both  ends  of  the  branch.  The  pulse 
generator  was  connected  to  one  wire  at  junction  J^.  The  results  to  be 
studied  will  be  the  transient  voltages  across  the  resistors  at  junction 


The  three  wires  have  different  radii  and  different  dielectric 
coatings,  as  shown  in  Figure  lb.  The  centers  of  the  wires  are  all 
4.16  cm  above  the  ground  plane,  and  the  total  length  of  the  transmission 
line  is  20  meters. 

The  transmission  line  parameters  were  obtained  by  using  the 
time  domain  reflectometry  measurement  technique  [Ref.  3] .  The  per-unit- 
length  capacitance  and  inductance  matrices  were  measured  to  be 


46.4807 
[C]  =  -20.9093 

-20.5534 


[L]  = 


0.883944 

0.484157 

0.535475 


-20.9093 
33.8339 
-  4.1526 


0.484157 

0.939843 

0.378936 


-20.5534 

-  4.1526  pF/m 
31.0993 

0.535475 
0.378936  yH/m 
0.992230 


The  characteristic  impedance  matrix  of  the  line  is  given  by 


232.419 

148.382 

159.830 


148.382 

237.879 

121.279 


159.830 

121.279 

258.320 
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Since  there  are  three  wires  coated  with  three  different  dielectric 
materials,  there  are  three  separate  quasi-TEM  modes  traveling  at 
different  velocities.  The  modal  velocities  are  2.829  *  10®  m/s, 

2.009  *  10®  m/s,  2.178  *  10®  m/s,  respectively. 

The  open  circuit  voltage  of  the  pulse  generator  was  measured 
and  is  shown  in  Figure  2.  The  pulse  has  a  peak  of  about  6  volts. 

The  half-power  pulse  width  is  about  2  ns. 

Table  2  lists  the  input  data  cards  for  QV7T  code. 

Figures  3a,  b,  c.  Figures  4a,  b,  c,  and  Figures  5a,  b,  c  show 
the  transient  responses  at  junction  Jg  when  one  wire  at  junction  Jj 
is  driven  by  the  pulse  generator.  Both  the  calculated  and  measured 
results  are  plotted  in  the  figures,  and  they  all  show  very  good 
agreement.  It  is  apparent  in  the  figures  that  there  are  three  trans¬ 
mission  line  modes  traveling  at  different  velocities  on  the  transmission 
line. 

It  should  be  noted  that  the  results  shown  here  do  not  include 
secondary  reflections  at  the  loads.  It  took  the  fastest  mode  70  ns  and 
the  slowest  mode  100  ns  to  travel  from  to  Jg.  The  first  secondary 
reflection  pulse  at  O2  has  t0  Travel  3  *  20  meters  and  it  occurs 
about  210  ns  after  the  turn-on  time  of  the  pulse  generator.  The 
measurements  were  not  taken  for  such  a  long  duration  of  the  response. 

The  secondary  reflection  pulses  are  usually  low  in  magnitude,  hence 
the  sensitivity  and  noise  cause  problems  in  the  measurement.  The 
calculations  were  not  performed  to  obtain  the  secondary  reflections 
because  it  is  necessary  to  reduce  the  time  resolution  in  order  tc 
accommodate  the  larger  time  window. 

For  most  of  the  responses  shown  in  the  comparisons,  it  was 
found  that  the  amplitudes  of  the  calculated  pulses  are  in  general 
larger  than  the  measured  values.  This  could  be  attributable  to  the 
transmission  line  losses  in  the  experiment. 
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TABLE  2.  LISTING  OF  INPUT  DATA  FOR  3-WIRE  TRANSMISSION  LINE  SECTION 


$TOPO  NB=1,NJ=2$ 

$BRNH  IB=1,NWW=3$ 

$JNC  JN=1,NBOO=1,IBB(1)=1$ 

$JNC  0N=2,NBJ0=1,IBB(1)=1$  %  „  „  „  0 

$LDPTS  NL0AD=3,LL0AD(1)=1»1,1»0L0AD(1)=2,2,2,LLN0(1)=1,2,3$ 
$LINEPRM  XK(1,1)=46.4807E-12,-20.9093E-12,-20.5534E-12, 

XK( 1 ,2 ) =-20 . 9093E- 12 ,33 . 8339E- 12 , -4 . 1526E- 12 , 

XL ( 1 . 1 ) =0/993944E-6 ,0 . 484157E-6 . 0 . 535475E-6 , 

XL( 1 ,2)=0 . 484157E-6 ,0 . 939843E-6 ,0 . 378936E-6 , 

XL( 1 ,3)=0 . 535475E-6 ,0 . 378936E-6 ,0 . 992230E-6 , 
XLEN=20 .0 , IFLAG=0$ 

$SOURCES  XLP0S=0.0,VSRC(1)=1. 0,0. 0,0.0  Note  1. 

$CCNN  JN=1,N0CC=1$ 

$LDC  JZL=2$ 

$LOAD  2(1)=25. 0,50.0, 50.0$  Note  2. 

$CONN  JN=2,N0CC=1$ 

$LDC  JZL=2$ 

LOAD  Z(l)=50. 0,50. 0,50.0$ 

$FFTDR  TMAX=400. ,M=10$ 

$WVFORM  ISTEP=3$  „  ,  „ 

$WVMATCH  NTT=13,TT( 1)=0.0,2.4,3.5,4.4,5.6,6.8,7 .8, 
9.0,11.0,12.4,14.4,16.0,18.0, 

VTT(1)=0. 0,3. 0,6. 0,3. 0,-0. 1,0. 3, 0.05, 0.25, 
0.1,0.2,0.0,0.18,0.0,SCAL=1.0 


Note  1.  If  wire  2  at  is  driven,  change  to 
VSRC(1)=0. 0,1. 0,0.0$ 

If  wire  3  at  J1  is  driven,  change  to 
VSRC(1)=0. 0,0. 0,1.0$ 

Note  2.  If  wire  2  at  is  driven,  change  to 
Z(l)=50. 0,25. 0,50.0$ 

If  wire  3  at  is  driven,  change  to 
Z(l)  =  50.0,50.0,25.0$ 
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Figure*  'll).  Trauu lent  voltage  of  wire  2  at  J„  when  wire  i  at  ta  driven. 
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figure  3c.  Transient  voltage  of  wire  3  at  J2  when  wire  1  at  Is  driven. 
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figure  4c.  Transient  voltage  of  wire  3  at  J„  when  wire  2  at  J.  is  driven 
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figure  5a.  Transient  voltage  of  wire  l  at  J„  when  wire  3  at  Jj  is  driven. 
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figure  lb.  Transient  voliagc  of  wire*  2  at  when  wire  3  at  .1  is  driven. 


B.  T-network 

The  T-network  being  investigated  involves  a  junction  of  multi¬ 
conductor  transmission  lines.  The  network  geometry  is  shown  in 
Figure  6a  and  the  wire  configurations  are  shown  in  Figures  6b,  c,  and  d. 

The  network  topology  involves  three  branches  denoted  by  B^, 

B2,  Bg,  and,  four  junctions  denoted  by  J^,  J2>  Jg,  J^.  These  three 
branches,  respectively,  have  five,  three  and  two  wires.  The  wires  at 
junction  were  loaded  with  50ft  resistors  to  the  ground.  Wire  1  at 
junction  1  was  driven  by  a  voltage  source.  The  wires  of  branch  2  were 
all  terminated  at  junction  Jg  with  50ft  resistors  to  the  ground.  The 
wires  of  branch  3  were  terminated  at  junction  J4  also  with  50ft 
resistors  to  the  ground.  The  three  branches  are  connected  at 
junction  J2  with  the  branching  of  wires  as  shown  in  Figure  6a. 

The  voltage  sources  considered  are  a  pulse  generator  and 
a  unit-step  function  generator.  The  open  circuit  voltages  of  the 
generator  are  shown  in  Figure  7.  The  internal  resistances  of  the 
generators  are  50ft  as  shown  in  the  network  geometry  of  Figure  6. 

The  transient  voltage  responses  of  the  resistors  at  junctions  Jg  and 
J4  have  been  measured  and  will  be  compared  with  calculated  results. 

The  configurations  of  the  wires  of  all  branches  are  shown 
in  Figures  6b,  c,  d.  The  radius  of  the  conductors  and  insulators  and 
the  dielectric  constants  of  the  insulators  of  various  wires  are  listed 
in  Table  3. 

The  per-unit-length  capacitance  and  inductance  matrices  were 
obtained  from  the  time  domain  reflectometry  measurements.  Thus,  the 
[C]  and  [L]  matrices  for  the  5-wire  line  were  measured  to  be. 


”52.555 

-16.125 

-17.761 

-8.895 

-5.632’ 

-16.125 

42.145 

-2.319 

-15.357 

-14.264 

= 

-17.761 

-2.319 

35.779 

-3.145 

-6.841 

-  8.895 

-15.357 

-3.145 

41.771 

-9.705 

_-  5.6325 

-14.264 

-6.841 

-9.705 

25.355 

m* 

I 
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Figures  6b,  c  and  d.  The  wire  configurations  of  the  T-network 


TABLE  3.  WIRE  RADIUS  AND  DIELECTRIC  CONSTANTS 


Conductor  radius 

Insulator  radius 

i 

*  0.108  cm 

»  0.35  cm 

i 

a2  •  0.062  cm 

b2  ■  0.33  cm 

a^  *  0.062  cm 

*  0.15  cm 

a.  *  0.111  cm 

4 

n 

a^  *  0.0125  cm 

*  0.25  cm 

^  -  2.25 


e2  -  4.10 


2.7  <  e3  <  3.0 


s4  -  1.65 


Eg  *  1.41 


0.902 

0.4775 

0.5399 

0.48516 

0.5133 

0.4775 

0.9513 

0.38426 

0.48609 

0.4239 

0.5399 

0.38426 

1.0067 

0.4247 

0.5343 

0.48516 

0.48609 

0.4247 

0.8507 

0.5568 

0.51334 

0.4239 

0.5343 

0.5568 

1.116 

The  corresponding  characteristic  impedance  matrix  of  the  5-wire  line 
is  given  by 


223.073 

138.964 

150.885  , 

138.278 

145 

.372 

138.964 

229.164 

115.579 

140.970 

125 

.528 

150.885 

115.579 

247.950 

123.761 

149 

.056 

138.278 

140.970 

123.761 

225.051 

155 

.647 

145.372 

125.528 

149.056 

155.647 

295 

•  409_ 

The  line  parameters  of  the  3-wire  line  branch  B2  are  exactly  the  same 
as  those  given  in  the  3-wire  line  section  discussed  before.  The  per- 


unit-length  capacitance  and  inductance  matrices  of  the  2-wire  line 
in  branch  are 


tC] 

[L] 


P 26.746 

-14.684 

[-14.684 

20.052. 

I"  0.8295 

0.5334 

L 0.5334 

1.1012 

pF/m 


yH/m 


The  corresponding  impedance  matrix  of  the  2-wire  line  in  branch  B^  is 
given  by 


tZ0] 


227.59 

156.939 


156.939 

302.948 


0 


The  input  data  for  the  QV7T  program  for  this  T-network  is  shown 
in  Table  4.  The  voltage  responses  at  the  load  points  at  and 
have  been  computed  for  the  purpose  of  comparing  with  the  experimentally 
observed  data. 
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TABLE  4.  LISTING  OF  INPUT  OATA  FOR  THE  T-NETWORK 


STOPO  N8*3,NJ*4S 
SBRNH  IB*I,NWW*5$ 

SBRNH  I8*2,NUW«3S 
SBRNH  IB*3,NWW*2S 
SJNC  JN-1,NBJJ»1,IBB(1)*1S 
SJNC  JN*2,NBJJ*3,IBB(1)*1,2,3S 
SJNC  JN-3,NBJJ*1,IBB(1)*2S 
SJNC  JN*4 , N8J J*1 , IBB ( 1 ) *3S 

SLDPTS  NLOAO-5.JLOAD(l)«3,3,3,4,4,LLOAD(l)*2,2,2,3,3,LLNO(l)*l,2,3.1.25 
SLINEPRM  XK(l,l)*52.555E-12,-16.125E-12,-17.761E-12,-8.895E-12,-5.6325E-12, 

XK( 1,2)*-16.125E-12,42.1453£-12,-2.319E-12,-15.357E-12,-1.4264£-12, 

XK(l,3)*-17.761E-12,-2.319E-12,35.779E-12,-3.145E-12,-6.8411E-12, 

XK(l,4)*-8.895E-12,-15.357E-l2,-3.145E-12,41.771E-l2,-9.705£-I2, 

XK(1,5)*-5.6325E-12,-1.4264E-12,-6.8411E-12,-9.705E-12,25.355E-12, 

XL(l,l)*0.902E-6,0.4775E-6,0.5399E-6,0.48516E-6,0.51334E-6, 

XL(1,2)*0 .4775E-6 ,0 .9513E-6 ,0 . 38426E-6 ,0 . 48609E-6 ,0 . 4239E-6 , 
XL(l,3)*0.5399E-6,0.38426E-6,1.0067E-6,0.4247E-5,0.5343E-6, 

XL ( 1 ,4)*0.48516E-6 ,0. 48609E-6 ,0. 4247E-6 ,0 . 8507E-6 ,0 . 5568E-6 , 
XL(l,5)«0.51334E-6,0.4239E-6,0.5343E-6,0.5568E-6,1.115E-6, 

XLEN«8.0,IFLAG«0S 
SSOURCES  XLPOS*0.0,VSRC(1)*1.0S 

SLINEPRM  XK(1,1)*46.4807E-12,-20.9093E-12,-20.5534E-12, 

XK(  1 ,2)-20 .9093E-12 ,33 .8339E-12,  -4 . 1526E-12, 
XK(1,3)»-20.5534E-12,-4.1526E-12,31.0993E-12, 
XL(l,l)*0.883944E-6,0.484157E-6,0.535475E-6, 
XL(1.2)*0.484157E-6,0.939843E-6,0.378936E-6. 

XL( 1 ,3 ) «0 . 535475E-6 .0 . 378936E-6 ,0 .992230E-6 , 

XLEN-12.0.IFLAG-OS 

SSOURCESJ 

SLINEPRM  XK(1,1)*26.746E-12,-14.684E-12. 

XK(1,2)-14.684E-12,20.052E-12, 

XL ( 1 , 1 ) -0 . 8295E-6 ,0 . S334E-6 , 

XL ( 1 ,2) *0 . 5334E-6 , 1 . 1012E-6 , 

XLEN«8.0,IFLAG«OS 

SSOURCESS 

SCONN  JN-1.N0CC-1S 
SLOC  JZL*2S 

SLOAO  Z(l)»50.0,50.0, 50.0, 50.0,50.0$ 

SCONN  JN*2,N0CC*5S 
SCNM  JC(1)*1,0,0,0,0, 1,0,0,0,05 
SCNM  JC(1)«0, 1,0,0,0,0,1,0,0,05 
SCNM  JC(1)»0,0,1,0,0,0,0, 1,0,0$ 

SCNM  JC(1)»0,0,0,1,0,0,0,0,1,0$ 

SCNM  JC(1)«0,0,0,0, 1,0, 0,0, 0,15 
CONN  JN*3,N0CC*1S 
SLOC  JZL*2S 

SLOAO  Z(1 >*50.0,50.0,50.0$ 

SCONN  JN«4,N0CC*1$ 

sloc  ja*2S 

SLOAO  Z(l)*50. 0,50.0$ 

SFFTOR  TMAX«200.,M*9$ 

SWVFORM  ISTEP-3S  > 

SWVMATCH  NTT*13,TT( 1)*0. 0,2. 4, 3. 5, 4. 4, 5. 6, 6. 8, 7. 8,  Note  that  for  Unit- 

9.0,11.0,12.4,14.4,16.0,18.0,  .step  response,  change  to: 

VTT(1)*0. 0,3. 0,6. 0,3. 0,-0. 1,0. 3,0. 05,0. 25,  SWVFORM  ISTEP-lS 
0.1,0. 2,0. 0,0. 18, 0.0,SCAL*1.25S  land  delete  SWVMATCH. ..S 

SWVFORM  ISTEP-OS 
SWVFORM  ISTEP*OS 
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Figures  8a,  b,  c,  d,  and  e  are  the  transient  voltage  responses 
of  a  pulse  input  driving  on  wire  1  at  Jj.  For  all  the  load  points,  the 
comparison  between  the  measured  and  the  calculated  results  show  good 
agreement.  For  the  response  on  wire  '4  at  J^,  the  measured  curve  has  a 

very  sharp  peak  at  about  75  ns.  The  duration  of  the  peak  is  about 

1  ns.  We  were  unable  to  obtain  this  peak  from  the  calculation  because 
the  program  limits  us  to  the  use  of  time  steps  no  finer  than  1  ns. 

We  were  later  informed  by  the  personnel  of  Mission  Research  Corporation* 
that  their  direct  time  domain  calculation  also  did  not  show  that  par¬ 
ticular  peak.  From  their  judgment,  the  peak  may  be  caused  by  the 
coupling  of  nearby  objects  during  the  measurements. 

The  unit  step  responses  of  the  T-network  has  also  been  compared. 
The  excitation  in  the  calculations  is  an  ideal  unit  step  function  driving 
wire  1  at  junction  J^.  The  excitation  of  the  laboratory  measurement  is 
shown  in  Figure  7b,  and  is  observed  to  have  a  rise  time  of  about  2  to 

3  ns.  Although  this  idealization  will  affect  only  the  early  time 

responses,  it  should  be  pointed  out  that  the  deviation  of  the  calcula¬ 
tions  from  the  experimental  data  is  not  significant  even  at  early  times. 

Figures  9a,  b,  c,  d,  and  e  show  the  comparison  of  calculated  and 
measured  transient  responses  on  the  load  points  at  junctions  J3  and  J^. 


*  A.K.  Agrawal ,  Mission  Research  Corporation,  personal  communication, 
June  1978. 


30 


70  75  80  85  90  95  100  105  110  115  120 

Time  (ns) 


70  75  80  35  90  95  100  105  110  115  120 

Tine  (ns) 


Figure  8b.  Transient  voltage  on  wire  2  at  J, 
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Figure  8c.  Transient  voltage  on  wire  3  at 
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Figure  8d.  Transient  voltage  on  wire  4  at 
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Figure  8e.  Transient  voltage  on  wire  5  at  J4 
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Figure  9c.  Unit  step  response  on  wire  3  at  J3 
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Figure  9e.  Unit  step  response  on  wire  5  at  J. 


C.  The  H-network 


Having  completed  the  comparison  of  numerical  calculations 
with  the  experimental  data  for  the  3-wire  line  section  and  the 
T-network,  we  now  turn  our  attention  to  a  somewhat  more  complicated 
transmission  network.  A  diamond  shaped  network  had  originally  been 
proposed;  but  owing  to  the  complications  in  the  laboratory  setup, 
this  was  later  changed  to  an  H-shaped  network,  in  consultation  with 
AFWl  personnel.  The  transient  voltage  responses  in  an  H-network 
involving  two  transmission  line  junctions  were  experimentally 
measured  and  the  results  provided  [Ref.  5]  for  the  purpose  of  comparison 
with  numerical  results.  The  H-network  as  shown  in  Figure  10a 
consists  of  five  single-wire  transmission  line  sections  of  1.524  m, 

3.048  m  and  4.572  m  in  length.  The  transit  time  of  the  signals  in 
each  section  is  about  5  ns,  10  ns,  and  15  ns,  respectively. 

The  terminals  of  the  network  are  loaded  by  50fi  resistors 
as  shown  in  Figure  10a.  The  wire  at  junction  is  driven  by  a 
unit  step  function  generator.  The  open  circuit  waveform  of  the 
function  generator  has  been  previously  shown  in  Figure  7b.  The 
voltages  at  junctions  J^,  J^,  Jg,  and  dg  were  experimentally 
observed  and  numerically  computed  as  well. 

All  the  transmission  line  sections  have  the  same  configura¬ 
tions  and  a  typical  cross  section  of  the  lines  is  shown  in  Figure  10b. 
The  wires  are  supported  over  a  ground  plane  by  suitably  placed 
styrofoam  (er  -1)  blocks.  The  ground  plane  is  0.915  m  wide 
which  will  be  approximated  to  be  an  infinite  ground  plane.  The 
per-uni t-length  capacitance  of  the  wire  is* 


The  input  data  for  the  QV7T  program  are  shown  in  Table  5. 


*  EMP  Interaction  Handbook,  Section  2.4.2,  To  be  published  by  the 
Air  Force  Weapons  Laboratory,  Kirtland  AFB,  NM. 
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TABLE  5.  LISTING  OF  INPUT  DATA  FOR  THE  H-NETWORKS 


$TOP0  NB=5,NJ*6$ 

$BRNH  IB-1, NWW-1S 
$BRNH  IB=2,NWW=1$ 

SBRNH  IB*3,NWW=1$ 

$BRNH  IB=4,NWW=1$ 

$BRNH  IB=5,NWW=1$ 

$JNC  JN=1,NBJU=1,IBB(1)=1$ 

$JNC  JN=2,NBJU*3,IBB(1)=1,2,3$ 

$JNC  JN*3,NBJJ=3, IBB(1)=3,4,5$ 

$JNC  JN=4,NBJ0=1,IBB(1  =2$ 

$JNC  JN=5,NBJU=1 , IBB( lj=4$ 

$JNC  JN=6,NBJJ=1,IBB(1)=5$ 

$LDPTS  NLOAD=4,JL0AD(1)=1, 4,5,6, 

LL0AD(1)=1 ,2,4,5,  LLN0(1)=1,1,1,15 
“SLINEPRM  XK(1,1)=2.2697E-11,  XLEN=3.048,  IFLAG=1$ 
SSOURCES  XLP0S=0.0,  VSRC(1)=1.0,  CSRC(1)=0.0$ 

-  $LINEPRM  XK(1,1)=2.2697E-11,  XLEN=1.524,  IFLAG=1$ 
$SOURCES$ 

$LINEPRM  XK(1 ,1)=2.2697E-11,  XLEN-3.048,  IFLAG=1$ 
$SOURCES$ 

$LINEPRM  XK(1,1)=2.2697E-11,  XLEN=3.048,  IFLAG=1$ 
$SOURCES$ 

$LINEPRM  XK(1,1)=2.2697E=11,  XLEN-4.572,  IFLAG=1$ 
$SOURCES$ 

$CONN  JN=1$ 

$LDC  JZL=1$ 

$LOAD  Z(l)=50.0$ 

$CONN  JN=2,NOCC=l$ 

$CNM  JC(1)=1,1,1$ 

$CONN  JN=3  NOCC=l$ 

$CNM  JC(1)-1,1,1$ 

$CONN  JN=4$ 

$LDC  JZL*1$ 

$LOAD  Z(l)=50.0$ 

$CONN  0N=5$ 

$LDC  JZL*1$ 

$LOAD  Z(l)=50.0$ 

SCONN  JN=6$ 

SLDC  JZL*1$ 

$LOAD  Z(l)*50.0$ 

$FFTDR  TMAX=400 ,M=9$ 

SWVFORM  I STEP* 1$ 

SWVFORM  ISTEP*0$ 

SWVFORM  ISTEP*0$ 

SWVFORM  ISTEP*0$ 

SWVFORM  ISTEP*0$ 

TOPO  NB*0$ 


Figures  11a,.  b,  c,  d  show  the  transient  voltage  responses 
across  the  5012  resistors  at  every  termination  of  the  network.  As 
seen  In  Figures  lib  and  d,  the  measured  signals  are  found  to  travel 
slower  than  the  speed  of  light.  The  slowdown  may  be  caused  by  the 
actual  laboratory  environment  such  as  the  styrofoam,  the  finite 
ground  plane  and  other  nearby  objects. 
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TIME  (ns)  ,  ,  ,  . 

Figure  11a.  The  unit  step  response  at  junction  Jj(voUage  across  the  resistor) 
of  the  ll-network 


025 


The  unit  step  response  at  junction  J4  of  the  H-network 


TIME  (nt) 

response  at  junction  J5  of  the  H-network 


0.30 


Summarizing  this  section,  we  find  that  there  is  good  agreement 
between  the  measured  transient  response  on  the  laboratory  models  and 
the  numerically  evaluated  responses  which  make  use  of  the  QV7T  com¬ 
puter  code.  Since  the  measured  driving  waveform  for  the  pulse  input 
was  digitized  for  use  in  the  numerical  computations,  the  agreement  between 
the  theory  and  experiment  is  seen  to  be  somewhat  better  than  for  the 
step  input,  for  which  an  idealized  waveform  was  used  in  the  calculation. 
However,  it  should  be  pointed  out  for  the  case  of  the  step  Input 
the  departure  in  calculated  responses  from  the  measured  responses  is 
not  significant. 

In  the  following  section,  we  turn  our  attention  to  compute 
and  compare  the  transient  response  on  a  power  cable  inside  the  TACAMO 
aircraft,  with  experimentally  measured  data  in  the  ATHAMAS  I  (HPD) 
facility  at  AFWL. 


SECTION  III 


COMPARISON  WITH  TACAMO  POWER  CABLE  MEASUREMENT 

In  addition  to  the  laboratory  models  discussed  previously,  we 
have  also  obtained  adequate  data  on  some  cables  in  an  aircraft  from 
the  TACAMO  Add-On  Testing.  Specifically,  the  60-90  kVA  power  cable 
assembly  13499  which  is  connected  to  No.  1  engine  will  be  considered 
in  this  report  for  data  comparisons. 

The  power  cable  consists  of  seven  wires  for  the  three  phase 
power  system.  It  runs  parallel  to  the  edge  of  the  wing  from  the  engine 
and  enters  the  aircraft  interior  through  two  receptacles  located  at  the 
pressure  break.  The  two  bundles  of  the  cable  inside  the  airplane  body 
are  terminated  in  the  distribution  box. 

During  the  EMP  simulation  test,  it  is  conjectured  that  the 
power  cable  picks  up  energy  mainly  from  the  portion  of  the  cable  along 
the  edge  of  the  wing.  The  excitation  sources  of  the  cable  in  the  wing 
are  distributed  along  the  line,  and  hence  difficult  to  determine.  In 
order  to  simplify  the  measurement  and  analysis  we  shall  use  the 
Thevenin  equivalent  circuit  to  represent  the  outside  section  of  the 
cable  by  input  impedances  and  open  circuit  voltages.  Neglecting  the 
coupling  of  the  energy  into  the  power  cable  in  the  interior  of  the 
aircraft  body  from  any  other  sources,  we  are  able  to  establish  a  network 
model  with  the  equivalent  discrete  sources  at  the  pressure  break. 

Several  major  modifications  on  the  QV7T  computer  code  have 
been  made  to  enable  the  calculation  of  the  TACAMO  power  cable  analysis. 
They  are: 

1.  The  original  code  is  suitable  only  for  networks  with 
frequency  independent  loads.  The  program  was  modified 
to  include  the  cases  when  some  of  the  loads  are 
frequency  dependent. 

2.  Subroutines  were  added  to  generate  the  impedance  of  the 
equivalent  circuit  from  the  measured  curves. 

3.  The  code  was  modified  to  accept  the  numerical  input  wave¬ 
form  for  each  wire.  The  data  were  stored  in  a  magnetic 
tape. 

The  network  topology,  impedance,  sources,  transmission  line  parameters, 
and  the  results  of  the  data  comparisons  will  be  discussed  separately. 
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A.  Network  Topology 

The  relative  location  of  the  power  cable  in  the  aircraft  is 
shown  in  Figure  12.  The  cable  consists  of  seven  wires  including  one 
neutral  and  six  others  of  the  three  phase  power  system.  Thus,  each 
phase  has  two  wires  that  are  electrically  in  parallel. 

Figure  13  sketches  the  geometry  of  the  cable  inside  the  air¬ 
craft.  The  power  cable  actually  includes  two  bundles,  one  with  three 
wires  and  another  with  four  wires  (including  the  neutral).  The 
neutral  wire  branches  from  the  four-wire  bundle  and  is  grounded  * 
midway  between  the  pressure  break  and  the  distribution  box. 

The  terminations  of  the  cable  are  very  complicated  impedance 
networks.  For  the  purpose  of  comparing  with  the  computer  model, 
transmission  lines  with  simple  terminations  are  desired.  Thus,  the 
terminal  at  the  distribution  box  is  shorted  and  grounded  to  avoid  the 
complications  of  the  impedance  inside  the  distribution  box.  The 
termination  at  the  pressure  break  will  be  represented  by  the  Thdvenin 
equivalent  circuit  which  includes  the  excitation  sources.  The  impedance 
and  driving  waveforms  of  the  equivalent  circuit  will  be  discussed  in  the 
following  subsections. 

The  electrical  configuration  of  the  network  becomes  the  one 
shown  in  Figure  14.  The  junctions  and  branches  are  defined  as  follows: 

Junction  J^:  At  the  pressure  break  where  it  connects  to  the 
Th^venin  circuit. 

Junction  J2:  At  the  point  where  the  neutral  wire  branches  out  from 
the  bundle.  It  is  about  2  meters  in  length  from  the 
pressure  break. 

Junction  J^:  At  the  distribution  box  (grounded).  The  wires  reaching 
this  junction  are  connected  to  the  original  four  wire 
branch.  (Only  three  wires  reach  this  junction  because 
the  neutral  wire  branches  out.) 

Junction  J^:  At  the  distribution  box  (grounded).  The  wires  reaching 
this  junction  come  from  the  original  three  wire  bundle. 
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Figure  12.  The  location  of  the  power  cable  to  be  tested 


Branch  1:  This  branch  is  2  meters  long  and  includes  all  the  seven 
wires  from  to  J2- 

Branch  2:  This  branch  is  1.3  meters  long  and  includes  three  wires 
from  J2  to  Jj- 

Branch  3:  This  branch  is  1.0  meters  long  and  includes  three  wires 
from  Jg  to 

Branch  4:  This  branch  is  1.3  meters  long  and  has  only  the  neutral 
wire  from  J2  to  J5  (ground). 


B.  Equivalent  Impedance  at  Pressure  Break  Toward  Engine  No.  1. 

The  Thevenin  equivalent  circuit  requires  the  measurement  of 
the  impedance  matrix  for  the  portion  of  the  cable  outside  the  pres¬ 
sure  break.  A  network  analyzer  (HP8407A)  with  the  calibrated 
impedance  probe  (HP1165A)  was  used  to  measure  the  impedance.  The 
measurements  were  done  in  the  frequency  range  of  1  -  110  MHz. 

Since  the  impedance  probe  (HP1165A)  could  not  reach  the 
connectors  at  the  pressure  break  during  the  impedance  measurement, 
a  coaxial  cable  was  used  to  connect  between  the  probe  and  the  terminal 
Thus,  the  measured  impedance  curves  must  account  for  the  section  of 
the  coaxial  cable  using  the  Smith  chart  or  transmission  line  analysis. 
The  characteristic  impedance  of  the  coax  is  known  to  be  50ft.  Let 
the  wavenumber  and  the  length  of  the  coaxial  cable  be  B  and  l 
respectively.  Assume  that  St  is  proportional  to  the  frequency  f, 
we  have  ,  Kf 

where  K  (rad/MHz)  is  a  constant  independent  of  frequency.  The 
value  of  K  can  be  obtained  from  the  curve  of  the  100ft  terminated 
coax  as  shown  in  Figure  15.  Adjusting  the  curve  to  the  calibration 
values  given  by  the  dashed  lines  in  Figure  15,  we  found  that 

(1)  at  f-  30.59  MHz,  BJL  -  \ 

Hence,  K  z  (|)/30.59  z  0.05135  rad/MHz 

(2)  at  f  =  62.80  MHz,  &SL  *  n 

Hence,  K  =  tt/62.80  =  0.05003  rad/MHz 

(3)  at  f  =  90.90  MHz,  Bit  - 

Hence,  K  a  (|*)/90.90  =  0.05184  rad/MHz 
Taking  the  average  value  of  the  above  cases,  we  have 

St  -  0.05107  x  f 
where  f  is  the  frequency  in  MHz. 


ibratlon  Curve  of  the  Network  Analyser 
ual  measurement  curve 
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All  the  impedance  measurements  should  be  adjusted  according 
to  the  coax  parameters.  Figure  16  illustrates  the  setup  of  the 
measurement  with  coax.  From  the  classical  transmission  line  theory, 
we  have 


Z 


a 


Zm  cos  BA  -  j  Zc  sin  BA 
1Q  cos  BA  -  j  Zm  sin  BA 


Hence  the  actual  impedances,  Z  ,  are  obtained  from  the  measured 
impedances,  Zm>  using  the  above  equation. 

The  impedance  matrix  of  the  Thlvenin  equivalent  circuit  for 
the  power  cable  at  the  pressure  break  has  a  dimension  of  7  x  7.  A 
complete  characterization  should  include  28  individual  impedance 
measurements,  because  the  matrix  is  symmetric.  Studying  the  geometry 
of  the  three-phase  power  cable  between  the  pressure  break  and  the  engine 
as  shown  in  Figure  17,  we  found  that  only  eight  different  measurements 
are  required  because  of  the  similarities  between  wires.  Wires  A^,  Bj, 
C1  have  similar  configurations,  as  do  wires  A£,  B2,  C2.  The  impedance 
relations  can  be  divided  into  the  following  eight  groups  with  each 
group  having  identical  impedances. 


Group  1: 

A1 

-  Ground, 

Bl- 

Ground , 

ci 

-  Ground 

Group  2: 

-  Ground 

Group  3: 

A2 

-  Ground, 

B2- 

Ground, 

C2 

-  Ground 

Group  4: 

A1 

'  Bl*  A1 

'  cr 

B1 

-C1 

Group  5: 

A2 

-  B2,  A 2 

C2* 

B2 

C2 

Group  6: 

A1 

■  °r  Bi 

-v 

C1 

-D1 

Group  7: 

A2 

■  °r  b2 

■  Dr 

C2 

-D1 

Group  8: 

A1 

“  ^2 

B2* 

A1 

-C2 

81 

-  A2.  B1 

B2’ 

B1 

C2 

C1 

-  A2.  Cj 

B2* 

C1 

-C2 

• 

I 
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Figure  16.  Impedance  measurement  using  a  coaxial  cable. 


pressure  break 


Figure  17.  Power  cable  geometry  outside  the  pressure  break. 


(actual 

impedance) 


A  total  of  10  redundant  measurements  has  justified  the  above  groupings. 

Figures  19  through  26  are  the  corrected  impedance  values  of 
the  equivalent  circuit.  The  data  are  all  digitized  and  used  as  input 
to  the  modified  QV7T  program. 

The  recovered  impedance  curve  of  the  100ft  load  is  shown  in 
Figure  18  for  the  purpose  of  checking  the  accuracy  of  the  impedance 
measurement.  It  has  been  found  that  the  measurements  of  the  impedance 
values  are  not  very  accurate  at  high  frequencies.  The  approximation 
on  ,  the  stray  capacitances  of  the  end  of  the  coax,  and  the  quality 
of  the  resistor  may  be  some  of  the  reasons  for  the  discrepancies. 

At  frequencies  above  65  MHz,  the  phase  of  the  impedance  becomes  negative 
while  for  a  pure  resistor  it  should  remain  zero.  The  amplitude  of 
the  impedance  is  also  larger  than  100ft  for  frequencies  higher  than 
65  MHz.  The  maximum  phase  error  is  about  50°  at  110  MHz  and  the 
maximum  impedance  amplitude  error  is  about  50ft  at  110  MHz.  Hence, 
using  the  impedance  data  for  frequencies  higher  than  65  MHz  should  be 
avoided,  if  possible. 
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Figure 


Figure  19.  Power  cable  impedance  at  pressure  seal  looking  toward  the  wing 
(for  A,  -  Cround,  B  -  Ground,  C,  -  Ground) 


impedance  at  pressure  seal  looking  toward  the 


Figure  21.  Power  cable  impedance  at  pressure  seal  looking  toward  the  wing 
(for  A  -  Ground,  B„  -  Ground,  and  C_  -  Ground) 


toward 


Power  cable  impedance  at  pressure  seal  looking  toward  the  wing 
(for  A.  -  B„,  A„  -  C„,  B  -  C  ) 


le  impedance  at  pressure  seal  looking  toward  the  wing 


Power  cable  Impedance  at  pressure  seal  looking  toward  the  wing 


C.  The  Open  Circuit  Voltages  of  the  Thevenin  Circuit 


The  open  circuit  voltages  of  the  cable  at  the  pressure  break 
under  the  EMP  testing  are  the  input  sources  for  the  Thevenin  equivalent 
circuit.  In  simulation  testing  it  is  usually  easier  to  measure  the 
currents  than  the  voltages.  The  open  circuit  voltage  can  be  obtained 
by  measuring  the  current  in  the  configuration  of  Figure  27  and  by  the 
following  discussions. 

While  all  other  wires  are  open,  the  i-th  wire  is  grounded  and 
the  current  Im  is  measured.  Let  the  impedance  matrix  be  denoted 

by  tZ.  .1.  The  open  circuit  voltage  V  .  is  simply 

i  jJ  01 


Voi  ~  ^i,i  *mi 

for  i  =  1,  2,  3,  ....  7.  Where  Z.  .  is  the  i-th  diagonal  element 

*  9  * 

of  the  impedance  matrix. 

The  currents  were  measured  using  the  DASET  system  which  sent 
the  signal  to  the  DASET  center  via  a  dielectric  waveguide.  The  current 
probes  were  the  clamp-on  type.  The  corrected  time  domain  data  were 
stored  in  two  magnetic  tapes.  We  have  obtained  the  tapes  from 
Rockwell  International  Corporation,*  and  have  converted  them  into 
tapes  which  are  compatible  with  the  Lawrence  Berkeley  Laboratory 
tape  readers. 

The  time  domain  current  data  are  fed  into  the  QV7T  code  and 
transformed  into  the  frequency  domain  data  by  FFT.  For  each  frequency 
point  the  current  data  are  multiplied  by  the  proper  impedance  values 
to  convert  the  voltage  sources  of  the  Thevenin  circuit.  Figures  23a 
through  34a  are  the  transient  current  measured  during  the  simulator 
testing.  The  frequency  spectra  of  these  currents  is  shown  in 
Figures  28b  through  34b. 

In  using  the  measured  time  domain  current  data,  we  have 
discovered  discrepancies  between  the  data  in  magnetic  tapes  and  the 
curves  presented  in  the  graphs.  The  time  scales  recorded  on  the  tapes 
are  0  to  1000  ns  while  the  graphs  for  the  same  measurements  have  the 
time  scales  of  0  to  4800  ns  and  0  to  4960  ns.  We  found  that  the  data 


*  fi.E.  Morgan,  Rockwell  International  Corporation,  personal  communica¬ 
tion,  January  1977. 
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on  the  magnetic  tapes  are  incorrect  because  their  fundamental  frequencies 
are  about  five  times  higher  than  that  of  the  actual  aircraft  resonant 
frequency  (~3.8MHz).  The  time  domain  current  data  were  obtained  from 
the  magnetic  tapes  by  correcting  the  time  scales. 


>00  1800  2400 

Time  (ns) 

transient  current  measurement 
pressure  break 


The  transient  current  measurement  on  wire  B,  of  the 
pressure  break 


•e  291 


30a.  The  transient  current  measurement  on  wire  C,  of  the  4-wire  bundle  at  the 
pressure  break  1 


(neutral)  of 
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The  transient  current  measurement  on  wire  AQ  of  the  3-wire  bundle  at  the 
pressure  break  c 
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Figure  34b.  The  short  circuit  current  spectrum  on  wire  C?  of  the  3-wire  bundle  at  the 
pressure  break  c 


D.  The  Transmission  Line  Parameters  of  the  Power  Cable 


The  seven  wire  branch  of  the  power  cable  consists  of  two 
twisted  bundles,  one  with  four  wires  and  the  other  with  three  wires. 
Since  the  two  bundles  are  not  packed  together,  we  shall  assume  no 
coupling  between  the  bundles.  The  cross  section  of  the  four  and 
three  wire  bundles  is  shown  in  Figure  35.  The  four  wire  bundle 
has  a  high  degree  of  symmetry  against  the  ground  plane;  thus  only 
the  configuration  in  Figure  35a  will  be  used  in  calculating  the  line 
parameters.  The  three  wire  bundle  has  two  entirely  different  con¬ 
figurations  against  the  ground  plane  as  shown  in  Figure  35  b  and  c. 
There  are  many  different  positions  of  the  wires  when  they  are 
twisting  along  the  path.  The  transitional  positions  are  disregarded 
because  the  average  values  of  the  per-unit-length  capacitance  and 
inductance  matrices  will  be  used. 

The  radius  of  the  conductors  of  each  individual  wire  is 
3.5  mm  and  the  radius  of  the  insulator  is  5.6  m.  The  relative 
dielectric  constant  is  assumed  to  be  2. 

A  separate  computer  code  using  the  multi -series  expansion 
method  [Ref. 6]  was  used  to  obtain  the  per-unit-length  capacitance 
matrix  of  each  wire  configuration. 

The  per-unit-length  capacitance  matrix  of  the  wires  of 
Figure  35a  was  found  to  be 


"78.508 

-26.905 

-19.563 

-25.872 

-26.905 

64.210 

-23.936 

-24.471 

-19.563 

-23.936 

124.120 

-18.583 

-25.872 

-24.471 

-18.583 

110.680 

The  elements  of  the  matrix  are  averaged  when  the  wires  are  twisted. 
The  new  per-unit-length  capacitance  matrix  of  the  twisted  four  wire 
bundle  becomes 


(b)  3-wire  bundle  (O  3-wire  bundle 

when  twisted 


Figure  35.  The  cross  section  of  the  power  cable; 

!:;%•;!  dielectric  medium  with  a  relative 
dielectric  constant  --  . 


» 
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"94.380 

-23.222 

-23.222 

-23.222" 

-23.222 

94.380 

-23.222 

-23.222 

-23.222 

-23.222 

94.380 

-23.222 

-23.222 

-23.222 

-23.222 

94.380 

pF/m 


The  per-unit-length  capacitance  matrix  of  the  wires  of 


Figure  35b  is 


110.67 

-18.611 

-25.990 

tCj]  = 

-18.611 

110.67 

-25.990 

-25.990 

-25.990 

62.383. 

and  that  for  Figure  35c  is 

'64.206 

-27.008 

-23.941 

tc2i  = 

-27.008 

64.206 

-23.941 

-23.941 

-23.941 

115.48 

pF/m 


pF/m 


The  elements  of  the  matrices  [C^]  and  [Cg]  are  averaged  to  obtain 
the  new  per- unit-length  capacitance  matrix  of  the  twisted  three  wire 


bundle  as  follows 


[C] 


87.936 

-24.247 

-24.247 


-24.247 

87.936 

-24.247 


-24.247' 

-24.247 

87.936 


pF/m 


The  inductance  of  the  transmission  lines  for  the  quasi -TEM 

modes  are  usually  independent  of  the  dielectric  constants  which 

insulate  the  wires.  Hence  we  can  assume  e  =  l  and  obtain  the 

r 

capacitance  matrices  as  described  above  for  this  new  £r .  The  per- 
unit-length  inductance  matrix  can  be  obtained  from  the  per-unit- 
length  capacitance  matrix  (e  =  1)  by  using  the  following  equation 

[L]  *  -y  where  c  is  the  speed  of  light, 

c 

The  resulting  per-unit-length  inductance  matrix  for  the  twisted 
four  wire  bundle  is 
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[LI 


0.267706 

0.088478 

0.088478 

0.088478 


0.099478 

0.267706 

0.088478 

0.088478 


0.088478 

0.088478 

0.267706 

0.088478 


0.088478 
0.088478 
0.088478 
0 . 267706_ 


pH/m 


and  that  for  the  twisted  three  wire  bundle  is 


[L] 


0.276717 

0.102229 

0.102229 


0.102229 

0.276727 

0.102229 


0.102229 

0.102229 

0.276717 


pH/m 


The  portion  of  the  neutral  wire  leading  to  the  ground  from 


junction  Jg  is  a  single  wire  and  ground  plane.  The  per-unit-length 
capacitance  and  inductance  of  this  single  wire  line  over  a  ground 


plane  are 

C  =  80.563  pF/m 
L  =  1.37918  nH/m. 
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E.  The  Results  -  Comparison  of  the  Calculated  and  the  Measured  Current 
at  the  Distribution  Box 

For  the  purpose  of  comparisons,  only  the  current  responses  on 
the  six  wires  of  the  three-phase  terminating  at  the  distribution  box 
will  be  considered.  The  original  network  connected  to  the  power  cable 
at  the  distribution  box  is  quite  complex.  In  order  to  simplify  the 
measurements  and  analytical  calculations,  the  wires  at  the  distribution 
box  were  shorted  to  ground. 

The  numerical  results  were  obtained  by  the  computer  model  using 
the  network  topology  described  in  section  A.  The  impedance  and  sources 
of  the  Thevenin  equivalent  circuits  were  measured  according  to  the  dis¬ 
cussions  in  sections  B  and  C.  The  per-unit-length  capacitance  and 
induction  matrices  of  the  transmission  lines  were  obtained  in  Section  D. 

Figures  36a/b,  37  a/b,  and  39a/b  through  41a/b  compare  the 
measured  and  calculated  current  spectra  of  phases  B^,  Ag.  B£  and  Cg. 
The  measured  and  calculated  transient  wire  currents  are  shown  in 
Figures  36c/d,  37c/d,  and  39c/d  through  41  c/d.  Some  inadequacies  have 
been  found  in  the  measured  results  as  stated  below. 

(1)  The  current  on  the  wire  of  phase  of  the  4-wire  bundle  was 
not  measured.  Thus,  only  the  calculated  curves  of  this  wire 
are  presented  in  Figure  38. 

(2)  The  measured  wire  currents  on  phase  A^  of  the  4-wire  bundle 
and  phase  A2  of  the  3-wire  bundle,  shown  in  Figure  36  and 
Figure  39,  are  apparently  incorrect.  There  are  two  reasons 
to  conclude  the  incorrectness:  (a)  The  main  resonances  at 
about  3.2  MHz  are  missing  for  both  cases.  All  other  wire 
current  measurements  on  the  power  cable  showed  a  main 
resonance  at  3.2  MHz;  (b)  The  time  window  on  the  transient 
response  of  phase  A^  is  only  up  to  432  ns  and  that  of 
phase  A ^  is  1760  ns.  The  wire  current  of  phase  A2  decayed 

to  negligible  values  at  about  700  ns.  All  other  measurements 
on  the  transient  wire  currents  of  the  power  cable  show  strong 
signals  even  after  1500  ns. 
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|l(jw)|  (amp/Hz)  |I(jw)|  (amp/Hz) 


Time  (ns) 

(c)  Measured  transient  current  (data  incorrect) 


Time  (ns) 


(d)  Calculated  transient  current 

Figure  36c, d.  Comparison  of  transient  wire  current  on  phase  A 1  of 
the  4-wire  bundle  at  the  distribution  box 


|I(jui)|  (amp/Hz)  IK  Jo)  |  (amp/Hz) 


|I(ju>)|  (amp/Hz) 


r 


Fi'iure  38a.  Wire  current  spectrum  on  phase  of  the  4-wire  bundle 
at  the  distribution  box  (Not  measured) 
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(c)  Measured  transient  current 
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(d)  Calculated  transient  current 

Figure  40c,d.  Comparison  of  transient  wire  current  on  phase  B-  of. 

the  3-wlre  bundle  at  the  distribution  box  c 
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(c)  Measured  transient  current 
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(d)  Calculated  transient  current 

Figure  41c ,d.  Comparison  of  transient  wire  current  on  phase  Cg  of 
the  3-wlre  bundle  at  the  distribution  box 


Hence,  only  the  measurement  on  phase  of  the  4-wire  bundle  and 
phases  Bg,  Cg  of  the  3-wire  bundle  can  provide  reasonable  comparisons 
with  the  numerical  results. 

The  measured  and  calculated  wire  currents  on  phase  B^  are 
shown  in  Figure  37.  In  the  frequency  spectrum,  both  the  measurement 
and  the  calculation  showed  a  main  resonance  at  about  3.2  MHz  with 
current  amplitude  of  the  order  of  6  *  10  A/Hz.  The  lowest  and 
highest  frequencies  we  have  used  in  the  FFT  calculations  are  about 
1  MHz  and  40  MHz  respectively.  The  measured  transient  response  was 
recorded  from  0  to  4800  ns.  The  calculations  were  performed  for  times 
up  to  3000  ns.  It  was  found  that  the  peaks  of  the  calculated 
transient  current  occur  at  about  the  same  time  as  the  measured  values. 
The  peak  values  calculated  are  about  twice  as  large  as  the  peaks  in 
the  measurement.  Since  we  only  used  3000  ns  as  the  time  window  for 
the  FFT,  the  time  domain  response  of  the  numerical  results  had 
experienced  the  aliasing  effect.  This  problem  can  be  solved  by 
increasing  the  time  window  in  the  FFT  process. 

The  comparisons  of  the  wire  currents  on  the  phases  and 
C2  of  the  3-wire  bundle  are  shown  in  Figure  40  and  Figure  41,  respec¬ 
tively.  Both  the  current  spectra  show  the  main  resonance  at  about 
3.2  MHz  and  the  amplitudes  are  very  close  to  the  measured  results. 

The  calculated  transient  currents  for  both  wires  Bj  and  Cg  are  shown 
from  0  to  1700  ns.  The  low  frequency  oscillations  and  the  peak 
values  of  the  numerical  results  for  both  wires  are  very  close  to  the 
measured  values.  The  calculated  peaks  are  a  little  higher  than  the 
measured  peak  values.  It  should  be  noted  that  the  responses  from 
1700  ns  to  3000  ns  were  not  plotted  because  the  aliasing  effects 
dominated  the  calculations  in  this  time  frame. 

Although  the  wire  currents  of  the  phases  B^,  Bg,  Cg  showed 
reasonable  agreement  between  the  calculated  and  experimental  results, 
the  agreement  is  not  as  good  as  we  found  in  the  comparison  of  the 
laboratory  models,  discussed  in  section  II.  Several  difficulties 
had  been  found  in  studying  the  data  comparisons  with  the  TACAM0 
power  cable  measurement,  and  they  are  listed  below. 
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(1)  The  pul se-to-pul se  variation  of  the  simulation  causes 
inconsistent  measured  data. 

(2)  The  measurement  of  the  Thevenin  equivalent  impedance  network 
was  not  very  accurate.  The  error  of  the  measurement  was 
illustrated  by  the  lOOfl  load  at  the  end  of  the  coax  as  shown 

in  Figure  18.  The  error  is  especially  large  when  the  frequency 
is  higher  than  about  65  MHz.  In  addition,  the  network  analyzer 
displays  nonlinear  effects  at  frequencies  lower  than  about 
4  MHz.  The  measurement  from  1  to  4  MHz  may  not  be  very 
accurate  while  the  main  resonant  frequency  is  3.2  MHz.  This 
might  contribute  a  certain  amount  of  error  in  the  calculations. 

(3)  It  was  proposed  originally  to  use  the  Norton  equivalent 
network  to  represent  the  power  cable  outside  the  pressure  break. 
The  short  circuit  currents  must  be  measured  when  all  the  wires 
are  grounded  simultaneously.  The  currents  were  mistakenly 
measured  by  grounding  only  one  wire  at  a  time.  This  leads  to 
the  Thevenin  equivalent  network.  The  open  circuit  voltage  of 
the  Thevenin  network  was  obtained  by  multiplying  the  impedances 
with  the  measured  current  spectrum.  The  errors  in  the  measured 
Impedances  can  also  cause  some  errors  in  the  driving  waveform 
represented  by  the  open  circuit  voltage. 

(4)  The  transmission  line  model  using  the  twisted  wires  over  a 
ground  plane  is  only  an  estimation  of  the  actual  configuration. 
The  per-unit-length  capacitance  and  inductance  matrices  obtained 
from  the  twisted  wire  model  may  deviate  from  the  actual  values 
because  of  the  approximations. 

(5)  Most  of  the  measured  transient  currents  decay  very  slowly 
within  3  us.  This  may  cause  difficulties  in  adjusting 
the  time  window  and  the  time  steps  In  the  FFT  algorithm. 

If  the  time  window  is  taken  larger  to  cover  all  the  oscillations, 
the  time  resolutions  will  be  affected  because  of  the  DIMENSION 
limitations  in  the  FFT  routines. 


101 


SECTION  IV 


CONCLUSIONS  AND  RECOMMENDATIONS 

The  analytical  and  experimental  transient  responses  of  example 
multi  conductor  networks  were  studied  in  this  report.  The  computer  code 
QV7T  which  utilizes  the  BLT  equation,  was  modified  and  used  in  perform¬ 
ing  the  numerical  calculations  for  the  analytical  results.  Measurements 
were  furnished  for  both  laboratory  models  of  multi  conductor  networks 
and  the  EMP  simulation  of  cables  on  an  aircraft. 

The  laboratory  transmission  line  network  models  involved  a  3-wire 
transmission  line  section,  a  T-network,  and  an  H-network.  Two  different 
driving  waveforms  were  considered,  i.e.,  the  pulse  and  the  unit  step 
functions.  It  was  found  that  the  analytical  results  agree  with  the 
experimental  measurements  very  closely  for  all  the  laboratory  network 
models.  These  comparisons  enhance  our  confidence  in  the  QV7T  code  for 
computing  transient  responses  of  certain  types  of  multi conductor  trans¬ 
mission  line  networks. 

The  measured  response  of  a  power  cable  in  the  TACAMO  aircraft 
during  simulator  testing  was  compared  with  analysis  of  the  EMP  transient 
response.  The  computer  results  show  the  same  resonance  frequency  as  the 
measured  response.  It  was  found  that  the  computer  analysis  tends  to 
over-estimate  the  peak  values  of  the  response.  Several  difficulties 
were  encountered  while  analyzing  the  power  cable  and  they  have  been  dis¬ 
cussed,  in  Section  III. 

Admittedly,  the  studies  of  this  report  are  still  incomplete  for 
a  complete  understanding  of  the  networks.  In  order  to  have  more 
comparisons  with  the  QV7T  code,  it  is  advisable  to  perform  additional 
laboratory  measurements  on  networks  which  include  frequency  dependent 
loads  and  which  include  T >rtps  in  their  topology.  It  is  also  suggested 
to  measure  some  canonical  Mlticonductor  transmission  lines  in  the 
simulator  environment  and  to  measure  the  transmission  line  responses 
under  distributed  source  excitations. 

For  an  extension  of  the  analytical  studies,  we  suggest  the  SEM 
study  of  simple  laboratory  transmission  line  models  using  the  BLT 
equation.  In  addition  it  would  be  useful  to  use  the  Prony  analysis  of 
experimentally  measured  responses  to  recover  SEM  frequencies  of  simple 
laboratory  transmission  line  models  and  compare  with  the  values  obtained 
from  the  BLT  equation  analysis. 
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